Inhomogeneity Induces Resonance Coherence Peaks in Superconducting 

Bi2Sr2CaCu208+5 
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In this paper we analyze, using scanning tunneling spectroscopy, the density of electronic states 
in nearly optimally doped Bi2Sr2CaCu208+i in zero field. Focusing on the superconducting gap, 
we find patches of what appear to be two different phases in a background of some average gap, 
one with a relatively small gap and sharp large coherence peaks and one characterized by a large 
gap with broad weak coherence peaks. We compare these spectra with calculations of the local 
density of states for a simple phenomenological model in which a 2^o x 2^o patch with an enhanced 
or supressed d-wave gap amplitude is embedded in a region with a uniform average d-wave gap. 

PACS numbers: 74.72.Hs, 74.50.+r, 74.25.-q 



One of the surprising features revealed by Scanning 
Tunneling Microscopy (STM) studies of the high Tc su- 
perconductor, Bi2Sr2CaCu208+5(BSCCO), is a pattern 
of patches of what appears to be two different phases, 
with significant differences in their electronic structures 
QiHilE^IEIll ' There are regions of relatively small local 
gap, A(r) ~ 25 — 35meV in which the peak in the local 
density of states (LDOS) at F = A is relatively sharp in 
energy and the peak height is very large. Other regions 
have a larger gap, with A(r) 50— 75meV and broad and 
small peaks (See Fig. P). It is tempting (as is widely as- 
sumed) to associate these very different electronic struc- 
tures with two different bulk electronic phases: the small 
gap regions, because they appear to have distinct co- 
herence peaks, are identified as regions of "good" su- 
perconductivity, whereas the large gap regions are like 
a pseudo-gap phase which competes with superconduc- 
tivity. This latter identification has found support from 
data suggesting that there is a subtle form of local charge- 
density wave order with period near four lattice con- 
stants ("stripes" or "checkerboards") 0, E i, S ESO 
which is most apparent in the large-gap regions [3, . 
However, because the characteristic size of the regions 
{L ~ 30A) is not much larger than the superconducting 
coherence length ( ^0 ~ 15A), it is clear that whatever 
the bulk character of each region, superconducting cor- 
relations can leak from one region into the other via the 
proximity effect j2|! thus complicating any such identifi- 
cation. 

In this paper we report results of STM studies on 
nearly optimally doped Bi2Sr2CaCu208+5 |i3 with high 
spatial and fine energy resolution. From this improved 
data we observe, as illustrated in Fig. that: 1) In 
the small gap regions, the peaks in the LDOS are too 
large to be the coherence peaks of a uniform ECS d-wave 



superconductor (see Fig. QJi) ; there is excessive spectral 
weight compared to the number of states pushed up from 
below the gap. 2) The peaks in the large gap regions are 
too broad and small to be the coherence peaks of a uni- 
form ECS d-wave superconductor (see Fig.^). 3) These 
regions are interspersed in a background "average" gap 
A = Aq w 40meV that produces a visible feature (typ- 
ically, a shoulder) in the LDOS in nearby regions; this 
coincides with the gap inferred from angle resolved pho- 
toemission measurements |l^ lisj . 

To interpret these results, we have calculated the quasi- 
particle LDOS for a mean field d-wave BCS model in 
which the strength of the pairing field (gap amplitude) is 
changed in a small L x L patch with L ~ 2^o- We find 
that structure like that seen in the "small gap regions" 
arises from resonant bound states if the gap amplitude 
vanishes (or is, at least, small compared to the peak en- 
ergy) in the Lx L patch, (see Fig.|3jD-) Structure similar 
to that seen in the large gap regions is found if a large 
pairing field is assumed inside such a patch jlS]. (see 
Fig. ISt.) One is thus led to conclude that small gap 
regions with "large coherence peaks" are regions with a 
much smaller than average pairing potential. Conversely, 
the fact that the concentration of large gap regions in- 
creases in increasingly underdoped samples, suggests that 
these regions, despite their strong pairing tendencies, 
have little or no superfluid density (phase stiffness) [l7j |. 
Finally, the fact that we see spectral features "leaking" 
between regions suggests that we are seeing patches of 
proximity coupled phases. 

For tunneling perpendicular to the Cu-0 planes, a typ- 
ical d-wave ECS shape of the spectrum is expected, char- 
acterized by a "v-shaped" LDOS at low bias and coher- 
ence peaks that accommodate the spectral weight from 
the opening of a gap with nodes. Early on, the general d- 
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wave shape of STM spectra was confirmed in BSSCO [l^ ; 
however data always appeared with significant particle- 
hole asymmetry in the background and subsequent anal- 
yses revealed that very few spectra quantitatively fit a 
BCS d-wave prediction, especially the coherence peak 
strength and shape. 

Figure shows spectra often seen in the small, aver- 
age, and large gap regions. In our analysis, we define 
the (positive bias) peak energy in the LDOS as the local 
gap A(r). In our samples, wc find the average gap is 
Aq 40 meV 0, 0, iS: iSl] and that approximately 75% of 
the area has a gap that is within ±10 meV of the average. 
Embedded in this background are patches of smaller gap 
(A < 30 meV) which cover approximately 15% of the 
area and of larger (A > 50 meV) which account for the 
remaining 10%. In addition to the earlier comments re- 
garding the coherence peak shape in these regions, we 
note that although the differences between the different 
spectra are more subtle 0, Q at energies far below the 
gap, the minimum tends to be more "v-shaped" in the 
large and average gap regions, and more rounded where 
the coherence peaks are anomalously large. 
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FIG. 1: Example spectra in a a) small, b) average, and c) 
large gap region. Inset is a typical 220A x 280A map of gap 
size. Crystal axes and squares that relate to Fig.|5|are marked. 

To study the behavior of the LDOS spectra as we go 
from a region of one gap size to another, we initially 
take a scan over a large area. Then we select several 
small areas and study them in detail, with a resolution of 
several spectra per atom. To maximize energy resolution, 
we limited the bias modulation used to acquire the dl/dV 
data to 2mV, and applied minimal filtering for the data 
collection. Figure El shows maps of A(r) with a small (a) 
and large (b) gap region. Below each figure we also show 
line cuts of spectra along the arrow. 

The spectra in Fig. [2^ illustrate the evolution of the 




FIG. 2: a) Gap distribution and the evolution of spectra as 
this region is crossed near a) small gap region, and b) large 
gap region. (Spectra follow arrows). 



LDOS on going from a small gap region with A < 30 
meV, to an "average" gap region with Aq ^ 40 meV. The 
crossover from one type of spectrum to the next occurs 
over a distance < with the anomalously large "coher- 
ence peaks" diminishing in strength while new peaks at 
the background gap of that region (typically ^ Aq) in- 
crease in intensity. The signature at Aq can be followed 
throughout the spectra from top (blue) to bottom (red): 
Deep in the small gap region it appears as a weak shoul- 
der above the main peak. At the border, a two-peak 
structure is apparent - one corresponding to the small 
gap characteristics and the other near Aq. Finally, out- 
side the small gap region, the main peak occurs at about 
Ao and has a more BCS-like structure. Line cuts through 
other small gap regions show similar behavior, although 
in instances of small gap differences, the two-peak struc- 
ture is harder to see. 

Figure |2b shows a line cut that starts in the average 
gap region, and ends up at the center of a large gap 
region. Again, the peak corresponding to the average 
gap diminishes in strength without dispersing in energy, 
while a new broad peak appears at A ~ 65 meV and 
the "average gap" coherence peak becomes a shoulder 
inside the large gap. Other line cuts in different large 
gap regions of the sample show similar evolution of the 
features, with the shoulder inside the gap being more/less 
visible when the differences in A are large/small. 

Finally, as can be noted from Fig. ^ (inset), while small 
gap regions may appear isolated in the background of 
the average gap, large gap regions almost always appear 
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within a distance < of small gap regions. This fea- 
ture may reflect the optimal-doping samples we are using 
which favor the creation of small gap regions. 

To obtain insight into the implications of these results, 
we have computed the LDOS for a simple model Hamil- 
tonian, meant to represent an effective mean-field Hamil- 
tonian for the quasiparticles: 



t.,s 



(1) 



where the vectors t label the lattice sites, 5 are the 
nearest-neighbor vectors, and {—Y = 1 for 5 — zti and 
{—Y = — 1 for 5 = ±y. In these expressions, we have 
adopted the usual Nambu notation with -0] = {c\^,cii). 
In the uniform case, A(i?) = Aq, 7i describes a uni- 
form square lattice with near-neighbor hopping t and 
a d-wave mean field characterized by a gap A(fc) — 
■^(cos fca; — cos ky). In the following we will consider the 
non-uniform situation in which A (^) = A on the sites in- 
side an L X L cluster embedded in a much larger (Af x M 
with M :^ L) cluster in which A(^) = Aq. All the cal- 
culations shown in the present paper are for M = 800 
and L — 5, although we have performed calculations for 
a range of Af 's and confirmed that M = 800 is large 
enough that the results are independent of M. 

We are interested in determining how the LDOS 
N{ll},£) varies as one moves from outside the cluster to 
sites inside the cluster, where 



N{uj,e):^ N{uj) (2) 
J2'Tr {G{e - e2)T{i2,h)G{h - £)) 



—Im 

TT 



Here N{lj) is the average density of states, runs over 
sites inside the L x L patch, G{£) is the single particle 
Green's function of the uniform lattice 



1 /cj + e^Ts + AfcTi 



(3) 



and T{£2, £i) is the T- matrix associated with the scatter- 
ing "potential" A{£) = A{£) ~ Aq, 



Ti£2,£i) = J2{-Y^nSi,,i,_s 
+E'(-)'^riG(£2 + 5-4m4,^i) . (4) 

Representative results of our calculations are shown in 
Fig. 3 for Ao/t = 0.2 and a small damping factor of 0.01. 



Results for the density of states at the center of a 5 x 5 
cluster in which the gap amplitude in the cluster ranges 
from A — 2Aq to A = are shown in Fig. 3a. Here, one 
sees that when the gap amplitude in the cluster is large 
compared to the background, the density of states at the 
center of the patch has a broad response at uj ~ ±A. 
However, as the cluster gap amplitude decreases below 
Ao, resonant peaks develop below ±Ao and move down 
in energy as A decreases. The height of the resonance 
peaks also increases as A decreases. In order to illus- 
trate the spatial dependence of N{uj,£), we consider the 
case in which A = 0, corresponding to a zero pairing am- 
plitude in the 5x5 patch. For this case, the density of 
states N{u}, £) versus uj for various sites [£x, iy) are shown 
in Fig. 3b. Here, the site (0, 0) corresponds to the center 
of the cluster and results are shown for £y = Q with 
varying from to 6. For sites inside the "gapless" cluster 
one sees a resonant response. This response appears at 
a lower energy than Aq and the peak height, for a given 
broadening, is significantly larger than the logarithmic 
structure in the bulk d-wave density of states. As one 
moves out from the center of the cluster, the sharp low 
energy peak in the LDOS loses intensity and weight be- 
gins to grow at w ^ Aq, so that at the boundary of the 
cluster a two-peak structure is clearly observed. Outside 
the cluster, the LDOS returns to its average behavior 
within a few lattice constants. 

If the gap parameter is doubled inside the cluster, 
A = 2Ao, one has iV(w, £) shown in Fig.|3t. In this case, 
for sites inside the cluster, N{uj,£) exhibits a broadened 
response near 2Ao as well as a weak response at Aq. Here 
again, as one moves several lattice spacings outside the 
cluster, the density of states returns to its uniform be- 
havior, characterized by the logarithmic coherence peaks 
at UJ — ±Ao. Note the change in scale between Figs.Otjb 
and 13; and how much stronger the resonance peaks are 
compared to the logarithmic peaks. 

The low energy behavior oi N{u}, £) is less dramatically 
i dependent than the peak structure. None-the-less, we 
believe that it is significant that the gap minimum is more 
"v-shaped" near the center of the large-gap cluster, and 
more rounded near the center of the small-gap cluster. 
We have not systematically explored the dependence of 
the results on the size of the cluster, L, but we have 
checked that similar behavior is obtained for somewhat 
different sized clusters and for clusters rotated by 45°. 

The model we have solved is admittedly overly simple, 
especially in that it neglects the strong on-site Hubbard 
U . Nonetheless, the qualitative similarities between fea- 
tures of the model calculations and the STM data suggest 
that some aspects of the problem are being successfully 
modeled. Clearly there are structural variations from 
place to place in Bi2Sr2CaCu208+5; for instance, a posi- 
tive correlation between the concentration of Oxygen in- 
terstitials and the large gap regions have been reported 
by McElroy et al. [lO|. However, can structural differ- 
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FIG. 3: LDOS N{uj, lx,iy) versus uj for a 5 x 5 patch centered 
at (0,0). a) LDOS at the center of the cluster N{uj, 0, 0) versus 
uj for different values of A(^). b) A(^) = on the patch, c) 
A(^) = 2Ao on the patch and for both, N{ijj,£x,£y) is shown 
for sites along {£x,0). 



ences due to these impurities give rise to a significant 
enhancement of the local pairing amplitude and what 
about the regimes in which the gap is suppressed? 

Rather, it would seem that some type of intrinsic am- 
plification of the effect of the structural variations is likely 
to be essential. For instance, if (as has been previously 
suggested doped antiferromagnets are near the cusp 
of a first order transition between two electronically dis- 
tinct states, then small differences in the local structure 
can nucleate small regions of one phase or the other. 
It has been well established by now that above opti- 
mal doping {x « 0.15 holes per Cu atom in the Cu-0 
plane) the average gap decreases with increasing dop- 
ing in proportion to T^, roughly as 2A « SksTc- In 
underdoped samples, the average gap increases with de- 
creasing doping, rising from Aq 40meV at optimal 
doping to Ao ^ 55meV at around x ~ 0.05, where 
Tc ^ Correspondingly, STM studies of under- 

doped Bi2Sr2CaCu208+5 reveal that the fraction of large 



gap regions increases with decreasing cc, and the fraction 
of small gap regions decreases 0,0,0. Thus, it would be 
natural to identify the large gap regions as being more 
representative of the electronic structure of underdoped 
cuprates and the small gap regions more representative 
of overdoped cuprates, both influenced by some "average 
gap" background. However, decreasing x also leads to 
a rapid decrease of Tc and the superfluid density, which 
implies that a large pairing field alone is insufficient to 
characterize the features of the electronic structure which 
refiect the approach to the Mott insulator. 
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